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By using our previous results of magnetohydrodynamical simulations for the solar wind from open flux tubes, 
I discuss how the solar wind in the past is different from the current solar wind. The simulations are performed 
in fixed one-dimensional super-radially open magnetic flux tubes by inputing various types of fluctuations from 
the photosphere, which automatically determines solar wind properties in a forward manner. The three important 
parameters which determine physical properties of the solar wind are surface fluctuation, magnetic field strengths, 
and the configuration of magnetic flux tubes. Adjusting these parameters to the sun at earlier times in a qualitative 
sense, I infer that the quasi-steady-state component of the solar wind in the past was denser and slightly slower if 
the effect of the magneto-centrifugal force is not significant. I also discuss effects of magneto-centrifugal force and 
roles of coronal mass ejections. 
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1. Introduction 

Young solar-type stars are generally very active: X-ray 
flux is up to ~ 1000 times larger than the present-day 
sun dGiidel et all Il997t iGudel L 12004. and the X-ray tem- 



peratu re is also higher (Rib as et all |2005; Tel leschi et al 
2005). Observations of young main sequence stars 
show very strong magn etic field strengths with an order 



of kG or even larger dDonati & Collier Cameron 1997 



Saar & Brandenburg 1 Il999[ ISaarl 12001 1 see also Donati & 



Landstreet 2009 for recent review), which are much stronger 
than the average strength of 1-10 G of the current sun. Young 
stars are generally fast rotators and strong magnetic fields are 
generated by strong differential rotat ion through dynam o ac- 



tivities of magnetoconvection (e.g., iBrun et all [2004). As 



time goes on, a star l oses its angular momentum through 
magnetic stellar wind ( Weber & Davis . 1967b a nd the mag- 
netic activities become weak ("e.g. JAvres Ul997l) . 

Solar wind is hot plasma emanating from the Sun, and the 
mass loss rate amounts to ~ 10 12 g s _1 (2 x 10 _14 M Q yr _1 ) 
at present. As inferred from high activity of young solar- 
type stars, I expect that the solar w ind was stronger at ear- 
lier times. IWood et al.l (|2002, 2005) observed asterospheres 
of low-mass stars by the Hubble space telescope and deter- 
mined the mass loss rates. The estimated mass loss rates 
show a decreasing trend with time, except for very young 
stars, thought there is still a large scatter in the data, partly 
because the observed stars range from very low-mass to near 
solar-mass starfl At the very early epoch, the mass loss rate 
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Although the observed stars include not only G-type (solar-analogue) stars 
but also lower mass K- and M-type stars, it is worth discussing trends be- 
cause the mechanism of mass loss of these stars are supposed to be similar; 
the stellar wind is accelerated by wave and turbulent motions associated 
with magnetic field driven by fluctuation in a surface convective zone. 



seems to be saturated at ~ 100 times of the present solar 
value, but some very young stars show lower mass loss rates 
than the present Sun. 

The solar wind roughly consists of two components. The 
first component is called fast solar wind with speed of 700 - 
800 km s _1 at the earth orbit. Fast solar wind is more steady 
and streams out from coronal holes which generally corre- 
sponds to open magnetic flux tube regions. The other com- 
ponent is slow solar wind with speed ~ 400 km s -1 . Slow 
wind is more complex and transient, and in most cases comes 
from lower latitude regions. Recent HINODE observations 
show that some portions of the slow wind appear to be origi 



nating from open flux tubes near active regio ns (ISakao et al 



|2007[ ilmada et all 120071: Itlarra et all 120081) . while there are 
still debates of other sources for slow wind. 

As other types of magnetized plasma from the sun, coro- 
nal mass ejections (CMEs) also are supposed to affect atmo- 
spheres of planets especially at early times. In this paper, 
however, I mainly focus on solar wind from open magnetic 
flux tubes, and briefly mention effects of CMEs afterward. 
Our group has carried out forward numerical modeling of the 
solar wind by using direct magnetohydrodynamical (MHD) 
simulations, in which we can directly test the response of the 
solar atmosphere to the surface fluctuations and properties of 
magnetic fields. By using the results of the simulations with 
parameters suitable for the younger sun, I discuss the solar 
wind evolution. 

2. Simulation Model 

In this section, I briefly describe the simulations of the 
solar wind. For detail, please refer to ISuzuki & Inutsuka 



(2006, SI06, hereafter). In order to cover the region with 
huge density contrast from p = 10~ 7 g cm~ 3 at the photo- 
sphere to p ~ 10~ 21 g cm~ 3 at the outer simulation bound- 
ary located at ps 0.1 astronomical unit (AU), simple one di- 
mensional (1-D) open flux tubes are adopted. The effects of 
super-radial expansion of flux tubes are incorporated by tak- 
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ing into account an expansion factor in the conservation of 
magnetic flux of radial (r) component, B r : 



B r r 2 f(r) — const., 



(1) 



where f(r) is a super-radial expansion factor. Most of the 
solar surface is covered by closed magnetic loop structure. 
Then, open flux tubes rapidly open above these loops; / 
is i ntroduced to consider this effect. The same function as 
in dKopp&Holzeiill976i) is adopted for f(r) (see SI06 for 
detail). 

In SI 06 we input the transverse fluctuations of the field 
line by the granulations at the photosphere. Amplitude, 
(dv±), at the photosphere is chosen to be compatible with 
the observed photosph eric velocity amplitude ~ 1km s _1 
dHolweger et al. . 1978h . SI06 tested various types of spec- 
tra; in this paper I discuss the results of the power spectra in 
proportion to 1/v, where v is frequency. The surface fluctua- 
tions generate upgoing Alfven waves which contribute to the 
acceleration of the solar wind in upper regions. At the outer 
boundaries, a non-reflecting condition is imposed for all the 
MHD waves, which enables us to carry out simulations for 
a long time until quasi-steady state solutions are obtained 
without unphysical wave reflection. 

SI06 dynamically treat the propagation and dissipation of 
the waves and the heating and acceleration of the plasma 
by solving ideal MHD equations with the relevant physical 
processes, including the sun's gravity, radiative cooling, and 
thermal conduction dSuzuki&Inutsukal 120051 SI06). I do 
not take into account stellar rotation, which may be impor- 
tant for very young stars. For the initial condition, we assume 
a static atmosphere with temperature T — 10 4 K in order to 
see whether the atmosphere is heated up to coronal tempera- 
ture and accelerated to accomplish the transonic flow. From 
t = 0, the transverse fluctuations are injected from the pho- 
tosphere and continue the simulations until the quasi-steady 
states are achieved. 

3. Solar Wind Response to Surface Condition 

By injecting fluctuations from the photospheric surface, 
the initially cool and static atmosphere is effectively heated 
and accelerated by the dissipation of the generated upgo- 
ing Alfven waves. The sharp transition region which di- 
vides the cool chromosphere with T ~ 10 4 K and the 
hot corona with T ~ 10 6 K is formed owing to a ther- 
mally unstable re gion around T ~ 10 5 K in the rad iative 
cooling function (Landi ni & Monsignori-Fossii Il990» . The 
hot corona streams out as the transonic solar wind (see 
Suzuki & InutsukaLl2005[ & SI06 for detail). 

The heating and acceleration of the solar wind plasma 
in inner heliosphere results from the dissipation of 
Alfven waves. In the simple 1-D treatment, Alfven waves 
mainly dissipate via nonlinear mode conversion; slow 
MHD waves are nonlinearly generated from outgoing 
Alfven waves and the slow MHD waves are damped by 
shock formation as a result of steepening of wave shape. The 
shocks also heat up surrounding plasma, which play a central 
role in the heating of the solar wind plasma. Magnetic pres- 
sure associated with the Alfven waves decreases with height 
as a result of the successive dissipation process, which di- 
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Fig. 1. Structures of corona and solar wind for different (dv±_) =0.4 
(blue), 0.7 (black dashed), and 1.4 km s _1 (red). I plot solar wind 
speed, v r (km s _1 ) (top left), temperature, T(K) (top right), density in 
logarithmic scale, log(p(g cm -3 )) (bottom right), and rms transverse 
velocity, (dvj_)(km s -1 ) (bottom left). Each variable is averaged with 
respect to time during 28min (the longest wave period considered). 



rectly pushes the plasma outward (momentum input) in ad- 
dition to gas pressure. 

The young sun is more active than the present-day sun, 
and therefore, I expect that surface fluctuations are stronger. 
The properties of the magnetic fields are also supposed to be 
very different. I discuss how properties of the solar wind are 
affected by fluctuation amplitudes and magnetic fields at the 
solar surface. 

3.1 Dependence on Surface Fluctuation Amplitude 

Figure Q] shows the response of the solar wind plasma to 
the surface fluctuations with different amplitudes, (dv± ) = 
1.4km s -1 , 0.7km s -1 and0.4kms _1 . The intermediate case 
((dv±) = 0.7km s _1 ) is a reference case which explains the 
present-day solar wind. 

The maximum temperature of (dv±) = 0.4km s _1 case is 
~ 5 x 10 5 K (TableQ]), which is cooler than the usual corona. 
The density is lower than the fiducial case by 1-2 orders of 
magnitude because the sufficient mass cannot be supplied to 
the corona by the chromospheric evaporation owing to the 
low temperature; the evaporation is drastically suppressed 
as T decreases, since the conductive flux (F c oc T 5 ^ 2 ^-) 
sensitively depends on T. As a result, the mass flux (pv r ) 
becomes much lower, by a factor of ^50, than that of the 
present solar wind. 

On the other hand, the larger (dv±) = 1.4km s _1 case 
gives larger coronal density, although the coronal tempera- 
ture is only slightly higher than the temperature of the ref- 
erence case. As a result of the dense corona, the mass flux 
of the solar wind is ~10 times larger (Table [TJ. This case 
shows complicated structure of temperature and density near 
the solar surface. The temperature starts to increase from a 
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(dv ± ) 




Wp.lAU 


W0.1AU 


(n p ti) iau 


1.4(km/s) 


1.2 x 10 6 (K) 


50(cm- 3 ) 


500(km/s) 


2.5 x lO^crrrV 1 ) 


0.7(km/s) 


1.0 x 10 6 (K) 


5(cm- 3 ) 


550(km/s) 


2.8 x 10 8 (cm- 2 s _1 ) 


0.4(km/s) 


0.5 x 10 6 (K) 


0.05(cirT 3 ) 


1200(km/s) 


6 x lO^crrrV 1 ) 



Table 1. Summary of model parameters and results. The first column shows the input amplitudes of the surface fluctuation. The second column presents 
the maximum temperatures in the simulation regions. The third column shows proton number density at 1 AU, which is calculated from n v at the outer 
boundary of the simulation by assuming the conservation of the mass flux (n p vr 2 f =const.) with constant speed from 0.1 AU to 1 AU. The fourth 
column shows the solar wind speeds at 0.1 AU (~ the outer boundary of the simulation region). The fifth column shows mass flux of the solar wind at 
1 AU, which is the product of the third and fourth columns. 



deeper location around r ~ 0.005i?o than the other cases. 
Thanks to this, the decrease of the density is slower (larger 
pressure scale height) so that the density around r = l.OLR© 
is two orders of magnitude larger than that of the reference 
case with (dv±) = 0.7kms _1 . However, the temperature de- 
creases slightly instead of showing a monotonical increase; 
it cannot go o ver the peak of the radiative coolin g function 
at T ~ 10 5 K dLandini & Monsignori-Fossil 1 1990t) because 
the radiative loss is efficient, owing to the large density. The 
temperature again increases from r ~ l.O3i?0 and above 
there the corona forms. The coronal density is larger than 
the that of the reference case as explained above. 

A striking feature is that the mass flux of solar winds 
sensitively depends on the input wave amplitude at the solar 
photosphere. Since the injected energy is proportional to the 
square of dv±, the difference of the input energy between 
the small amplitude case ((dv±) — 0.4km s _1 ) and the large 
amplitude case ((dv±) = 1.4km s _1 ) is ~ 12. However, the 
output mass flux of the large amplitude case is ~ 400 times 
larger than the mass flux of the small amplitude case. 

The sensitive behavior on the surface amplitude can 
be understood by reflection and nonlinear dissipation of 
Alfven waves (SI06). Because the heating is smaller, the 
temperature is lower in the smaller (dv±) cases. Then, 
the scale height becomes smaller and the density decreases 
rapidly. The Alfven speed, va{= B r / \fAirp), changes more 
rapidly due to the rapid decrease of density and the wave 
shape is largely deformed, which enhances the reflection. 
When the input wave energy decreases, a positive feedback 
operates; a smaller fraction of the energy can reach the coro- 
nal region. As a result, the density and mass flux of the solar 
wind becomes much smaller than the decreasing factor of the 
input energy. 

In addition to the effect of the wave reflection, the non- 
linear dissipation of the Alfven waves also plays a role in 
the sensitive dependence on the input wave energy (SI06). 
When the input wave energy becomes smaller, the density 
becomes smaller as explained above. Then, the nonlinear- 
ity of Alfven wave, dv±/vA, decreases, not only because the 
amplitude, dv±, is smaller but also because the Alfven speed, 
ua(oc 1/ \fp), is larger. Therefore, the Alfven waves do not 
dissipate and the heating is reduced, which further decreases 
the density; this is another type of positive feedback, which 
also results in the sensitive dependence of the density on the 
input wave energy. 
3.2 Magnetic Fields 

At earlier epochs the sun is expected to have possessed 
stronger magnetic field than today, as inferred from obser- 



vations of young stars dDonati & Collier Cameron , 1997 
Saar & Brandenburg 1 1 19991) The configuration of magnetic 
fields as well as field strength plays an essential role in 
determining the physical solar wind conditions. In the 1- 
D treatment as adopted in this paper, a super-radial ex- 
pansion factor, /^determines the geometry of flux tubes. 
Wang & Sheeievl dl990i fl99lT) showed that the solar wind 
speed at ~ 1AU is anti-correlated with / from their long- 
term observations. This can be naively understood by ener- 
getics consideration; since / directly determines an adiabatic 
loss of plasma in flux tubes, solar wind plasma in largely ex- 
panding (larger /) flux tubes undergoes larger energy loss. 

By the comparison of the outflow speed obtained by in- 
terplanetary scintillatio n measurements wit h observed pho- 
tospheric field strength iKoiima et al.1 (120051) have found that 
the solar wind velocity is better correlated with surface mag- 
netic field strength, divided by the expansion factor, B r ^J f; 
They claim that B r ^/ f is a better index for the speed of so- 
lar wind than individual 1// or B r fl. The correlation be- 
tween Brxt/ f and solar wind speed i s well explain ed by 



Alfven waves in expanding flux tubes (ISuzuki , 2006). Be 



cause the energy flux of input Alfven waves is proportional to 
B ri o, the positive correlation of solar wind speed with B r $ is 
quite natural. The negative correlation with / reflects the adi- 
abatic loss of Alfven waves in expanding tubes as explained 
above. 

One can understood the correlation more specifically from 
nonlinear dissipation of Alfven waves as well. The non- 
linearity of the Alfven waves, (dv±)/vA is controlled by 
va oc B r oc B r> o/f in the outer region where the flux 
tube is already super-radially open. Wave energy does not 
effectively dissipate in the larger B r ^j f case in the sub- 
sonic region because of relatively small nonlinearity, and 
more energy remains in the supersonic region. In general, 
energy and momentum inputs in the supersonic region re- 
sult in higher wind speed, while those in the subsonic region 
raise the mass flux (pv r ) of the wind by increasing density 
(ILamers & CassinelliL 1 19991) . This indicates that the solar 
wind speed is positively correlated with B r $ / f. 

4. Evolution of Solar Wind 

I have described based on our numerical simulations that 
the solar wind from open magnetic flux tubes is mainly 
controlled by the three parameters: surface fluctuation am- 
plitude, (dv±), radial magnetic field strength at the photo- 
sphere, B r xi, and super-radial expansion factor, /, of flux 
tubes. Thus, if I can determine these parameters for the 
young sun, I can estimate the properties of solar wind at that 
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time. 

4.1 Surface Fluctuation Amplitude 

I can speculate that (dv±) was larger at earlier times from 
circumstantial evidences. First, the rotation of the sun should 
be faster at earlier times because of loss of angular mo- 
mentum with time. Hence, in addition to surface convec- 
tion, interior differential rotation will lead to stronger surface 
fluctuations. From observational facts, younger stars shows 
higher X-ray activities as discussed in the introduction sec- 
tion, which also anticipates larger (dv±). An important issue 
is that a small increase of (dv±) leads to a large increase of 
mass flux of solar wind. 

4.2 Magnetic Field -B rfi Szf- 

B r fi is also supposed to be larger than to- 
day as implied from o b servat io ns of young stars 



is slower, wh i ch is consistent with the observational trend 
(IKoiima et all l2005). because of larger adiabatic loss in open 



1997; Saar & Brandenburg . 



magnetic flux tubes. 

Based on these considerations, namely (i) the magnetic 
field structure of the young sun more resembles the current 
solar maximum than the solar minimum, and (ii) open flux 
tubes show larger expansion during the solar maximum be- 
cause larger surface area is occupied by closed structure, I 
infer that open magnetic flux tubes of the young sun have 
much larger /. I speculate that, even though the magnetic 
field strength at the footpoints is stronger, B r ,Q / f, is smaller 
at earlier time. 

4.3 Solar Wind in the Past 

I would like to discuss how the properties of solar wind 
change in the past based on the simulation results. I have 
stated that (dv±) was larger and B r ^/ f was smaller at ear- 
lier times. These conditions imply considerably denser but 
slightly slower solar wind in the past, although it is quite dif- 
ficult to quantitatively determine the physical condition. 

4.4 Limitations 
In this paper I neglect the effect of ma g neto-c e ntrifugal 



dDonati & Collier Cameron 

19991) . Rossby number, Rq, which is defined as the ratio 
of stellar rotation period to convective turn-over time, is an 
important index wh ich describes stellar magnetic activities 
( Noves et all 0.984). Surface magnetic field strength, 
-Bph, multiplied by a magnetic areal filling factor, F, is 

well-correlated with Rq 1 . Stars with small Rq < 0.1, which force in driving winds (IWeber & Davis L 119671 iMestel 
corr espond to fas t rotating young stars, have B^F = 1-10 
kG dSaarlbOOlh . In some stars .F's are obtained and the 
values are 50-70%, which shows that B p ^ is an order of 10 
kG in these stars. 

If all the magnetic flux is open to the interplanetary space, 
I can set / = 1 /J 7 ; in reality, however, closed loop structure 
is more dominant, and then, it is expected that / ^> 1 jT. 
I speculate that super-radial expansion, /, of the young sun 
was larger as discussed from now, although it is not simple 
to pin down the specific value. Recently, the configuration of 
magnetic field has been obse rved in a number of low-mass 
stars (e.g. lDonati et al. , 2008 ). Fast rotating solar mass stars 
possess non-axisymmetric poloidal co mponents with sub 



stanti al toroidal fields (see Fig. 3 of iDonati & Landstreet 



2009). This implies that the field configuration is quite com 
plicated, and is very different from ordered dipole structure. 
In comparison with 1 1 -year periodicities of present-day solar 
activity 0, the properties of magnetic field of the earlier sun 
resemble the condition of the present-day solar maximum 
rather than the solar minimum. 

During the solar minimum periods, the solar wind con- 
sists of fast wind which is from mid- to high-latitude regions, 
and slow wind from lower-latitude regions. On the other 
hand, during the solar maximum periods, most of the re- 



1 1968b . Under the typical solar condition, if the rotation pe 
riod is 4 days (6-7 times faster than the present sun) or less, 
magneto-centrifugal force plays an imp ortant role and th e 
wind speed becomes significantly higher dNewkirkl. fl 980). 

Recently. iHolzwarth & Jardine I (12007) calculated the evo- 
lution of solar wind by considering megneto-centrifugal 
force. They concluded that the solar wind in the past was 
faster but not so much denser than today, wh i ch are different 
from my conclusion. IHolzwarth & Jard ine (2007) adopted 
the density at a 'reference point' ,1.1 times of the stellar ra- 
dius, in a parameterized way; they assume power-law de- 
pendence on stellar rotation rate. Since the basal density 
str ongly constrains the mass loss rate (see, e.g. section 3 
of lLamers & Cassinelli . 1999) and the adopted dependence 
in their standard case is rather weak, the derived mass loss 
rate shows weak dependence with time. In a strict sense, the 
density at the coronal base, accordingly the density at the 
reference level, is determined by the energy balance among 
heating, thermal conduction, and radiative cooling (SI06; see 
also Suzuki 2007 for different types of stars). The base den- 
sities of young active stars are larger than those of quiescent 
stars, and the mass loss rates are supposed to be significantly 
larger. In this case the wind speed becomes slower because 
more mass needs to be pushed away. The sa me tendency on 



gions are occupied by transient slow wind (McC omas et al.L the co ronal base density is also reported by IStereborg et al 



2008). At the solar maximum magnetic fields are dominated 
by a dipole component, while at the solar minimum higher 
multi-pole moments become more important showing more 



complex field configuration (Hak amada et aU 120051) . Be 



cause the surface is mostly covered with closed magnetic 
fields during the solar maximum, the field strength in the 
outer regions is not as strong as that at the solar minimum 
although the field strength at the photosphere is stronger. 
This fact shows that / is larger during the solar maximum. 
The speed of the solar wind during the solar maximum 



(2011) who performed 3D MHD simulations for a young 
sun. 

It is not still possible to discuss the wind speed in a quan- 
titative sense at the moment. My conclusion of the slightly 
slower wind in the past might be corrected because magneto- 
centrifugal force could be significant. On the other hand, I 
suppose that the wind speed is not so high as estimated by 



Holzwarth & Jardine (2007) bacause the coronal base den- 



sity should larger at that time than their assumption. Accord- 
ingly, the de nsity of the solar wind could be larger than those 



calculated in Holzwarth & Jardine (2007) 



I would like to note that the 11 -year duration is modulated with time d e- 
pending activity even in last several hundred years (Miyahara et al., 2008) 



In this paper I focus on quasi-steady-state solar wind from 
open magnetic flux regions, and do not consider CMEs. 
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While the quasi-steady-state component dominantly con- 
tributes to the total mass loss from the present-day sun than 
CMEs, at early times CMEs might play a more important 
role because of stronger magnetic activities. From the fos- 
sil record on the lunar surface, it i s inferred that the speed 
of solar wind wa s faster in the past dRav & Hevmann , 1980; 
Newkirk , 1980). I infer that the fossil record comprises ef- 
fects of both quasi-steady state component and CMEs. If 
strong CMEs dominantly contribute to the mass loss, the 
bulk speed of plasma flow from the sun would be recognized 
as higher. I speculate that the fossil inprint on the moon is 
connected with high CME activities at early epochs. 

5. Summary 

Based on our previous results of MHD simulations, I dis- 
cuss the evolution of the quasi-steady-state component of so- 
lar wind from open flux tubes. The properties of the solar 
wind are determined by the three parameters, surface fluc- 
tuation amplitude, magnetic field strength, and expansion of 
open flux tube. Referring to observational data of young low- 
mass stars, I suppose that the surface amplitude was larger 
and the ratio of field strength to flux tube expansion was 
smaller. Following this speculation, the solar wind in the 
past was dense and slightly slower than today, while if the 
effect of magneto-centrifugal force is taken into accout, the 
speed of the solar wind might be higher at very early epoch 
than my conclusion. CMEs are also supposed play a more 
important role in the past than today because of higher mag- 
netic activities. In order to understand the effects of solar 
magnetized plasma on the atmospheres of planets, we need 
to understand the evolution of CMEs as well as the evoltuion 
of quasi-steady-state solar wind. 
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